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Abstract 

The aim of this paper is to represent any continuous local martingale 
as an almost sure limit of a nested sequence of simple, symmetric random 
walks, time changed by a discrete quadratic variation process. One basis 
of this is a similar construction of Brownian motion. The other major 
tool is a representation of continuous local martingales given by Dambis, 
Dubins and Schwarz (DDS) in terms of Brownian motion time-changed 
by the quadratic variation. Rates of convergence (which are conjectured 
to be nearly optimal in the given setting) are also supplied. A neces- 
sary and sufficient condition for the independence of the random walks 
and the discrete time changes or, equivalently, for the independence of 
the DDS Brownian motion and the quadratic variation is proved to be 
the symmetry of increments of the martingale given the past, which is a 
reformulation of an earlier result by Ocone [8]. 

1 Introduction 

The present authors are convinced that both for theoretical and practical rea- 
sons, it is useful to search for strong (i.e. pathwise, almost sure) approximations 
of stochastic processes by simple random walks (RWs). The prototype of such 
efforts was the construction of Brownian motion (BM) as an almost sure limit 
of simple RW paths, given by Frank Knight in 1962 [7J. Later this construction 
was simplified and somewhat improved by Pal Revesz [S] and then by one of 
the authors It is interesting that this method is asymptotically equivalent 
to Skorohod embedding of a nested sequence of RWs into BM, [ITJ Theorem 
4]. This elementary approach also led to a strong approximation of Ito inte- 
grals of smooth functions of BM jTTJ Theorem 6]. We mention that this RW 
construction was extended to fractional BMs as well [T2] . 
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This paper discusses a generalization to continuous local martingales M. Be- 
side the above-mentioned RW construction of BM, the other major tool applied 
by the paper is a representation of continuous local martingales by Brownian 
motion, time-changed by the quadratic variation, given by Dambis pQ and Du- 
bins - Schwarz [3] (DDS). Theorem [2] shows that the quadratic variation process 
(M, M) can be almost surely uniformly approximated by a discrete quadratic 
variation processes N m which are based on stopping times of a Skorohod-type 
embedding of nested simple RWs into M. This corresponds to an earlier similar 
result by Karandikar [2] . Theorems [3] and 2] give an approximation of M by a 
nested sequence of RWs B m , time-changed by (M,M) and N m , respectively. 
The approximations almost surely uniformly converge on bounded intervals. 
Rates of convergence (which conjectured to be nearly optimal in the given set- 
ting) are also supplied. 

It is important to note that the DDS Brownian motion W and the quadratic 
variation (A/, M) are not independent in general, just like the approximating 
RW B m and the discrete quadratic variation N m . Since this could be a hindrance 
both in the theory and applications, a necessary and sufficient condition is given 
for the independence in Theorem [5] Namely, the approximating RW B m and 
the discrete quadratic variation N m (and so W and (M, M) ) are independent if 
and only if M has symmetric increments given the past. This is a reformulation 
of an earlier result by Ocone [8J, see [5] and [12] as well. 

Possible applications of the discrete approximations discussed in this paper 
include (a) generating continuous local martingales for which one has a suitable 
method to generate the discrete quadratic variation process N m , especially when 
the independence mentioned in the previous paragraph holds and (b) giving an 
alternative definition of stochastic integration with almost surely converging 
sums. (In a forthcoming paper the present authors are going to take up the 
second topic, extending some results of jllj.) 



2 Random walks and the Wiener process 

A main tool of the present paper is an elementary construction of the Wiener 
process (= BM). The specific construction we are going to use in the sequel, 
taken from [TT] , is based on a nested sequence of simple random walks that uni- 
formly converges to the Wiener process on bounded intervals with probability 1. 
This will be called RW construction in the sequel. One of our intentions in this 
paper is to extend the underlying "twist and shrink" algorithm to continuous 
local martingales. 

We summarize the major steps of the RW construction here, see [12] as 
well. We start with an infinite matrix of i.i.d. random variables X m {k) 1 
P {X m {k) = ±1} = 1/2 (m > 0, k > 1), defined on the same underlying prob- 
ability space (fi, T, P). Each row of this matrix is a basis of an approximation 
of the Wiener process with a dyadic step size At = 2~ 2m in time and a corre- 
sponding step size Ax = 2~ m in space, illustrated by the next table. 

The second step of the construction is twisting. From the independent ran- 
dom walks we want to create dependent ones so that after shrinking temporal 
and spatial step sizes, each consecutive RW becomes a refinement of the previ- 
ous one. Since the spatial unit will be halved at each consecutive row, we define 
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Table 1: The starting setting for the RW construction of BM 



At 


Ax 


i.i.d. sequence 


RW 


1 


1 


X (l),Xo(2),Xo(3) ) ... 


So(n) 


— 2^fc= 


-i^o(fc) 


2- 2 


2" 1 


X 1 (1),X 1 (2),X 1 (3) ) ... 


Si(n) 




=1*1 W 


2- 4 


2- 2 


X 2 (1),X 2 (2),X 2 (3),... 


5 2 (n) 




=1*2 W 



stopping times by T m (0) = 0, and for fc > 0, 

T m (k + 1) = min{n : n > T m (k), \S m (n) - S m {T m (k))\ =2} (m > 1) 

These are the random time instants when a RW visits even integers, different 
from the previous one. After shrinking the spatial unit by half, a suitable modi- 
fication of this RW will visit the same integers in the same order as the previous 
RW. We operate here on each point lu G ft of the sample space separately, i.e. 
we fix a sample path of each RW. We define twisted RWs S m recursively for 
k = 1, 2, . . . using S m -i, starting with So(n) = So(n) (n > 0). With each fixed 
m we proceed for k — 0, 1, 2, . . . successively, and for every n in the correspond- 
ing bridge, T m (k) < n < T m (k + 1). Any bridge is flipped if its sign differs from 
the desired: 

r n) = [ X m (n) if S m (T m (k + 1)) - S m (T m (k)) = 2X m ^{k + 1), 
[ — X m (n) otherwise, 

and then S m {n) — S m (n — 1) + X m (n). Then S m (n) (n > 0) is still a simple 
symmetric random walk [TTJ Lemma 1]. The twisted RWs have the desired 
refinement property: 

^S m {T m (k)) = Sm-tik) (m>l,fc>0). 

The last step of the RW construction is shrinking. The sample paths of 
S m {n) {n > 0) can be extended to continuous functions by linear interpolation, 
this way one gets S m (t) (t > 0) for real t. Then we define the mth approximating 
RW by 

B m (t)=2- m S m (t2 2m ). 
Using the definition of T m and B m we also get the general refinement property 

B m+1 (r m+1 (fc)2- 2 (™ +1 ') = B m (k2- 2m ) (m > 0,k > 0). (1) 

Note that a refinement takes the same dyadic values in the same order as the 
previous shrunken walk, but there is a time lag in general: 

T m+ i(fc)2- 2 ( m+1 ) - fc2~ 2m ^ 0. (2) 

Then we quote some important facts from QT| about the above RW con- 
struction that will be used in the sequel. These will be stated in somewhat 
stronger forms but can be read easily from the proofs in the cited reference, cf. 
Lemmas 2-4 and Theorem 3 there. 
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Lemma A. Suppose that X\, X2, ■ ■ ■ , Xn is an i.i.d. sequence of random vari- 
ables, E(Xfe) = 0, Var(Xfe) = 1. and their moment generating function is finite 
in a neighborhood of 0. Let Sj = X\ + • • • + Xj, 1 < j ' < N. Then for any 
C > 1 and N > N (C) one has 



sup \Sj\ > (2CN\ogN)i } < 2N X ~ C . 
i<j<N J 

We mention that this basic fact, that appears in the above-mentioned refer- 
ence [TT], essentially depends on a large deviation theorem. 

We have a more convenient result in a special case of Hoeffding's inequality, 
cf. 0]. Let Xl,X2, . . . be a sequence of bounded i.i.d. random variables, such 
that hi < Xi < en, and let S n — Y^=i-^i- Then by Hoeffding's inequality, for 
any x > we have 



/ n 

|S„-E(S n )|>x [j^di-h) 



4 ■ 
t=i 



< 2e~~. 



If E(Xi) = and b t = -a, here, then \ - hf = Y*=i A = Var(S„) if 

and only if X, = a t X[, where P {X[ = ±1} = ±, 1 < i < n. 

Thus if S 1 = a r X' r , where not all a r are zero and Var(S') = J^r a r < °°' 
we get 

P{\S\ > x(Vax{S))H < 2e~^ (x > 0). (3) 

The summation above may extend either to finitely many or to countably many 
terms. Let Si, S2, ■ ■ • , Sn be arbitrary sums of the above type: <Sfc = a krX' kr , 
P {X' kr = ±1} = i, 1 < /c < N, where X' kr and X' ls can be dependent when 
k 7^ Z. Then by the inequality ((3]) we obtain the following analog of Lemma lAl 
for any C > 1 and N > 1, 

sup |5 fc | > (2ClogiV)5 sup (Var(S fc ))' 

l<fc<JV l<fc<AT 
iV 

< ^Pjl^l > (2ClogiVVar( > 5'A ; ))^} < 27Ve- clogAr = 2iV 1 - c . (4) 
k=l 

Lemma |^] easily implies that the time lags ([2]) are uniformly small if m is 
large enough. 

Lemma B. For any K > 0, C > 1, and for any m > too(C), we have 
P < sup |T„ l+1 (fc)2- 2 (" i+1 ) - k2- 2m \ > ( lcK\og*K ) ' mh~ m I 

[o<A:2- 2m <X \2 / J 

< 2(X2 2m ) 1 ~ c , 

where log*x = max{l, logs}. 

This lemma and the refinement property ([T]) implies the uniform closeness 
of two consecutive approximations if m is large enough. 
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Lemma C. For any K > 0, C > 1, and for any m > m\{C), we have 
P< sup \B m+1 (k2- 2m ) -B m {k2- 2m )\ > Kf(log*K)im2- ! ? 

yo<k2- 2m <K 

< i{K2 2m f- 

where = max{l,if}. 

Based on this lemma, it is not difficult to show the following convergence 
result. 

Theorem A. The shrunken RWs B m (t) (t > 0, m = 0, 1, 2, . . . ) almost surely 
uniformly converge to a Wiener process W(t) (t > 0) on any compact interval 
[0, K], K > 0. For any K > 0, C > 3/2, and for any m > m^iC), we have 

P { sup \W{t) -B m (t)\ > Kl(log„K)im2- 2 ? 1 < G^ 2 " 1 ) 1 ' . 

[o<t<K J 

Now taking C = 3 in Theorem [3] and using the Borel-Cantelli lemma, we 

get 

sup \W(t) - B m (t)\ < 0(l)m2~^ a.s. (to -> oo) 

0<*<i<" 

and 

sup |W(t) -B m (t)| < Ki(\ogK)i a.s. (X ^ oo) 

0<t<i<" 

for any to large enough, to > TO2(3). 

Next we are going to study the properties of another nested sequence of ran- 
dom walks, obtained by Skorohod embedding. This sequence is not identical, 
though asymptotically equivalent to the above RW construction, cf. [TTJ The- 
orem 4]. Given a Wiener process W, first we define the stopping times which 
yield the Skorohod embedded process B m (k2~ 2m ) into W. For every m > let 
s m (0) = and 

Sm (k + l) = mf{s:s> Sm (k),\W( S )-W(s m (k))\ = 2- m } (k > 0). (5) 

With these stopping times the embedded process by definition is 

B m (k2- 2m ) = W(s m (k)) (to > 0, k > 0). (6) 

This definition of B m can be extended to any real t > by pathwise linear in- 
terpolation. The next lemma describes some useful facts about the relationship 
between B m and B m . These follow from [TTJ Lemmas 5,7 and Theorem 4], with 
some minor modifications. 

In general, roughly saying, B m is more useful when someone wants to gen- 
erate stochastic processes from scratch, while B m is more advantageous when 
someone needs a discrete approximation of given processes, like in the case of 
stochastic integration. 

Lemma D. For any C > 3/2, K > 0, take the following subset of the sample 
space: 

A m = \ sup sup |2- 2 ™T m ,„(fc) - k2- 2m \ < 6(CK* log, K)*mH- m \ , 

yn>m 0<k2- 2m <K J 

(7) 



G 



B. Szekely <fe T. Szabados 



where T m ^ n (k) = T n o T„_i o • • • o T m (k) for n > m > and A; > 0. TTien /or 
an?/ m > m^{C), 

P{A c m } < 4(X2 2m ) 1 " c . 

Moreover, lim n _»oo 2^ 2n T m>n {k) = t m {k) exists almost surely and on the set A m 
we have 

B m (k2- 2m ) = W(t m (k)) (0 < k2- 2m < K), 
cf. (OJ). Further, on A m except for a zero probability subset, s m (k) = t m {k) and 



sup \s m (k) -k2- 2m \ < 6{CK* log, K)imi 2~ 

0<k2~ 2m <K 



(m>m 3 (C)). (8) 



If the Wiener process is built by the RW construction described above using 
a sequence B m (m > 0) of nested RWs and then one constructs the Skorohod 
embedded RWs B rn (m > 0), it is natural to ask what the approximating 
properties of the latter are. The answer described by the next theorem is that 
they are essentially the same as the ones of B m , cf. Theorem Rl 

Theorem 1. For every K > 0, C > 3/2 and m > m^lC) we have 



P <j sup \W{t) - B m {t)\ > If* 4 (log* K)*m2-^ } < 10(K2 
n<t<K 



2m\l-C 



Proof. By the triangle inequality, 

sup \W(t) - B m (t)\ < sup W(t) - B m (t) 

0<t<K 0<t<K 



sup 

0<t<if 



B m (t) — B m (t) 



By Lemma [D] and equation (|6]) , on the set A m defined by Q we have 

B m (k2~ 2m ) = W{s m (k)) = B m (k2~ 2m ), 

except for a zero probability subset when m > 7713(C) . Since both B m (t) and 
B m (t) are obtained by pathwise linear interpolation based on the vertices at 
k2~ 2m E [0, if], they are identical on A m , except for a zero probability subset 
of it when m > m^C). Thus 



P sup \W(t)-B m (t)\ >KZ{\og*K)*m2-* 

{0<t<K 



< P{4 C „} + P<^ sup W(t)-B m (t) >K}(\og,K)*m2- 

lO<t<K 

Then by Theorem [X\ and Lemma ID] we get the statement of the theorem. □ 



3 The basic approximation results 

Beside the RW construction of standard Brownian motion, the other main tool 
applied in this paper is a theorem of Dambis (1965) and Dubins-Schwarz (1965) 
and an extension of it, cf. Theorems IBI and [Cl below. Briefly saying, these theo- 
rems state that any continuous local martingale (M(t),t > 0) can be trans- 
formed into a standard Brownian motion by time-change. Then somewhat 
loosely speaking, the resulting Brownian motion takes on the same values in 
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the same order as M(t), only the corresponding time instants may differ. These 
and other necessary matters about continuous local martingales will be taken 
from and discussed in the style of [TU] in the sequel. 

Below it is supposed that an increasing family of sub-cr- algebras {F t ,t > 
0) is given in the probability space {tl,F, P) and the given continuous local 
martingale M is adapted to it. 

In the case of a continuous local martingale M(t) vanishing at its quadratic 
variation (M, M) t is a process with almost surely continuous and non-decreasing 
sample paths vanishing at 0. This will be one of the two time-changes we are 
going to use in the sequel. The other one is a quasi-inverse of the quadratic 
variation: 

T s = inf {t : (M, M) t > s}, (9) 

where inf(0) = oo by definition. Then the sample paths of the process T s are 
almost surely increasing, but only right-continuous, since such a path has a jump 
at any value where the quadratic variation has a constant level-stretch. Beside 
this, T s may be infinite valued. The duality between the two time-changes is 
expressed by (M, M) t = inf {s : T s > t}. Observe that T s cannot have constant 
level-stretches since this would imply jumps for (M,M) t . Also the continuity 
of (M,M) t gives that (M,M) T , = s (s > 0), while we have only T {M ,M) t > t 
(t > 0) in the opposite direction. It is clear that 

(M, M) t < s=>t< T s , but t<T s =^> (M, M) t < s, (10) 

while 

(M, M) t <, >, > s t <, >, > T s , (11) 

respectively. 

Theorem B. [10, V (1.6), p. 181] If M is a continuous (Tt)-local martingale 
vanishing at and such that {MjM)^ = oo a.s., then W(s) — M(T S ) is an 
{^T s )-Brownian motion and M(t) — W((M,M) t ). 

Similar statement is true when (M, M)^ < oo is possible. Note that on the 
set {(M, M)oo < oo} the limit M(oo) = lim^oo M(t) exists with probability 1, 
cf. pig IV (1.26), p. 131]. 

Theorem C. [10, V (1.7), p. 182] If M is a continuous {IF t) -local martingale 
vanishing at and such that (M, M)^ < oo with positive probability, then there 
exists an enlargement (fi,.Ft,P) of {CI, J-y t ,P) and a Wiener process (3 on Q 
independent of M such that the process 

wn = ( _ ^a<(M,Af)oo, 

[S) \ M(oo)+,i3(s-(M,A/) 00 ) if s> {M, M). x 

is a standard Brownian motion and M{t) — W{{M, M) t ) for t > 0. 

From now on, W will always refer to the Wiener process obtained from M 
by the above time-change, the so-called DDS Wiener process ( or DDS Brownian 
motion) of M. 

Now Skorohod-type stopping times can be defined for M, similarly as for W 
in ©. For m > 0, let r m (0) = and 

T m {k + l)=mi{t:t>T m {k),\M{t)-M{T m {k))\=2- m } {k > 0). (12) 
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The (m+l)st stopping time sequence is a refinement of the mth in the sense that 
(•7"m(fc))felo * s a subsequence of (t„ 1+ i(j))°^ so that for any k > there exist 
ji and j 2 , T m+ i(ji) = T m (k) and r m+ i(j 2 ) = r m (k + 1), where the difference 
h - h > 2, even. 

Lemma 1. With the stopping times defined by from a continuous local 
martingale M one can directly obtain the sequence of shrunken RWs that almost 
surely converges to the DDS Wiener process W of M , cf. (£f|): 

B m (k2- 2m ) = W(s m (k)) = M(r m (fe)), s m (k) = (M, M) Tm(k) 

[but T m (k) < T Sm (u) ], where for m > 0, the non-negative integer k is taking 
values (depending on ui) until s m {k) < (M, M)^. 

Proof. By Theorems [B] and O it follows that W((M, M} Tm(fe) ) = M(r m (fc)). 
This implies that s m {k) < {M, M) Tm ^ k y Then consider first the case k = 1. If 
s m (l) < (M, M} Tm (!) held, then T Sm ^ < r m (l) would follow by dTTJ) , and this 
would lead to a contradiction because M(T Sm ( 1 )) = W(s m (l)) — ±2~ m . For 
values k > 1, induction with a similar argument can show the statement of the 
lemma. □ 

Below B m will always denote the sequence of shrunken RWs defined by 
Lemma [TJ 

Our next objective is to show that the quadratic variation of M can be 
obtained as an almost sure limit of a point process related to the above stopping 
times that we will call a discrete quadratic variation process: 

N m (t) = 2- 2m #{r:r>0 1 r m (r)<t} (13) 
= 2- 2m #{r :r >0,s m (r) < (M,M) t } {t > 0). 

Clearly, the paths of N m are non-decreasing pure jump functions, the jumping 
times being exactly the stopping times r m (fc). Moreover, N m (r m (fc)) = k2~ 2m 
and the magnitudes of jumps are constant 2~ 2m when m is fixed. 

Lemma 2. Let M be a continuous local martingale vanishing at 0, let (M, M) 
be the quadratic variation, T be its quasi-inverse |f5|). and N m be the discrete 
quadratic variation defined in \1S\) . Fix K > and take a sequence a m = 
0(m~ 2 ~ e 2 2m )K with some e > 0, where a m > K V 1 for any m > 1 (x V y = 
max(x, y), x A y = min(a:, y) ). 

(a) Then for any C > 3/2 and m > (C) we have 

sup \{M,M) t Aa m — N m (t A T am )\ > \2{Ca m log, a m )*m'2- 

0<t<K 

< Sia^ 2 ™) 1 - '. 

(b) Suppose that the quadratic variation satisfies the following tail- condition: 
a sequence (a m ) fulfilling the above assumptions can be chosen so that 

V{{M,M) t >a m }<D{t)m^~\ (14) 

where D(t) is some finite valued function of t € K+. Then for any C > 3/2 and 
m > m&(C) it follows that 

sup \{M,M) t - N m (t)\ > 12(Ca m log„a m ^m52- m 

0<t<K 

< 3{a m 2 2m ) l - c + D(K)m- x - 
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Proof. The basic idea of the proof is that the Skorohod stopping times of a 
Wiener process are asymptotically uniformly distributed as shown by ([5]) , while 
the case of a continuous local martingale can be reduced to the former by the 
DDS representation, cf. Lemma [1] 

i 1 

Introduce the abbreviation h a ^ m — 11.1 (Ca m log* a m ) 2 2~ m . Then h a , m 
= 0(to™ £ ) —> as to — > co. We need a truncation here using the sequence 
a m , since the quadratic variation (M,M) t is not a bounded random variable in 
general. By (IT) and ([T3|). 

JV m (*AT a J - 2- 2m #{r:r >0,r m (r) <tAT Qm } 

= 2~ 2m #{r : r > 0, s ro (r) < (M, M) t A a m }. 

On the event 

A a , m = I sup \s m {r) ~ r2- 2m \ < h a , m \ , 

[0<r2- 2m <2a m J 

if r = L((M, M) t Aa m +/i a ,m)2 2m J+l, then s m (r) > (M, M) t Aa m , so s m {r) is not 
included in N m {tf\T am ). Observe here that a m + h ajn + 2~ 2m < 2a m if m is large 
enough, m > (C), where we also suppose that 7714(C) > to 3 (C) and 7713(C) 
is defined by Lemma [Dj This explains why the sup is taken for r2~ 2m < 2a m in 
the definition of A a ^ m . Similarly on A a , m , if r — l({M,M) t A a m — h a , m )2 2m \, 
then s m (r) < (M, M) t A a m , so s m (r) must be included in N m (t A T am ). Hence 

(M,M} t Aa m -h a , m -2- 2m <N m (tAT a J < (M, M) t Aa m + h a , m + 2" 2m , (15) 

for any t G [0, if] on A a>m . 

Now 6(C2a m log >t (2 a „ l ))^TO5 2- m < 11.1 (Ca m log, a m )* m^2~ m = /i a , m , 
since log„(2a m ) < (1 + log2) log, a m . Hence it follows by Lemma [Dl that 

P {A c a , m } < 4(2a m 2 2m ) 1 - c < 3(a m 2 2m ) 1 - c , 

1 1 

when C > 3/2 and m > to 4 (C). Noticing that 0.9 (Ca m log, a TO ) 2 m^2~ m 
> 2~ 2m for any to > 1, this and (|T5"j) prove (a). 
Part (b) follows from (a), the inequality 

\(M,M) t -N m (t)\ < \{M,M) t -{M,M) t Aa m \ 

+ \(M, M) t A a m - N rn (t A T a J\ 

+ \N m (tAT a J-N m (t)\, (16) 

and from the following simple relationships between events: 

{(M, M) t A a rn + (M, M)t] = {(M, M) t > a m } 

and 

{N m (t AT a J ? N m (t)} = {t> T a J = {(M,M) t > a m }, 
cf. (dU). □ 
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Remark 1. When the quadratic variation {M,M) t is almost surely bounded 
above by a finite valued function g(t) > max{i, e} for each t>0, statements (a) 
and (b ) of Lemma [5| simplify as 

p( sup \(M,M) t -N m (t)\ >12(Cg(K)logg(K))imh- 

{0<t<K 

< 3(g(X)2 2m ) 1 - c 

for any K > 0, C > 3/2 and m > m^(C). 

The statement of the next theorem corresponds to the main result in Karand- 
ikar [5 , though the method applied is different and here we give a rate of 
convergence as well. 

Theorem 2. Using the same notations as in Lemma\^ and taking a sequence 
(c m ) increasing to oo arbitrary slowly, we have 

sup \(M,M) t - N m (t)\ < c m m^2~ m a.s. (m ->• oo). 

0<t<K 

Under the condition of Remark^ we also have 

sup \(M,M) t -N m (t)\ <g(K)i {log g(K))i a.s. {K -> oo) 

0<t<K 

for any m large enough, m > m^i). 

Proof. To show the first statement take e.g. C = 3/2 and a m = c rn in Lemma [5] 
(a). Consider the inequality (|T6|) . Since {M,M)k is finite-valued and c m — > oo, 
if to is large enough, depending on w, (M, M)k < a m holds and then t < T am 
holds as well by pU|) . These remarks show that the first and the third terms 
on the right hand side of inequality (|16[) are zero if to is large enough. Further, 
statement of Lemma [2] (a) can be applied to the second term. This, with the 
Borcl-Cantelli lemma, proves the first statement of the theorem. 

The second statement of theorem follows similarly from Lemma [5] (a) by the 
Borel-Cantelli lemma, taking C = 3 and a m = g(K). □ 

Now we are ready to discuss the strong approximation of continuous local 
martingales by time-changed random walks. 

Lemma 3. Let M be a continuous local martingale vanishing at 0, let (M, M) 
be the quadratic variation and T be its quasi-inverse (0). Denote by B rn the 
sequence of shrunken RWs embedded into M by Lemma\j^ Fix K > and take 
a sequence a m = 0(m~ 7 ~ e 2 2m )K with some e > 0, where a m > K V 1 for any 
to > 1. 

(a) Then for any C > 3/2 and to > m^(C) we have 
P I sup \M(t A T a J - B m ((M,M) t A a m )\ > fl&(tog„ o TO )*m2 _ * 

{o<t<K 

< 10(a m 2 2 ™) 1 - 

(b) Under the tail- condition J_?^[ ), for any C > 3/2 and m > m^iC) it follows 
that 

pj sup \M(t)-B m ((M,M) t )\>ai(\og,a m )im2-^ 

lO<t<K 

< 10(a m 2 2m ) 1 - c + Di^m- 1 -'. 
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Proof. First, take the DDS Wiener process W(s) obtained from M(t) by the 
time-change T s as described by Theorems [B] and [C] Since below we are going to 
use W(s) and also the time change T s only for arguments s < (M, M)^, we can 
always assume that W(s) = M(T S ) and M(t) = W((M,M) t ), irrespective of 
the fact whether (M, M)^ = oc or not. Second, define the nested sequence of 
shrunken RWs B m (s) by LemmaJTJ Then a quasi-inverse time-change (M, M) t 
is applied to B m (s) that gives B m ({M, M) t ) which will be the sequence of time- 
changed shrunken RWs approximating M(i). 
Since T s may have jumps, we get that 

sup \M(t)-B m ((M,M) t )\ > sup \M(T s )-B rn ((M,M) Ts )\ 

0<t<K Q<s<{M,M) K 

sup \W(s)-B m (s)\. (17) 

0<s<{M,M)k 

Recalling however that the intervals of constancy are the same for M{t) and for 
(M, M) t QI3 IV (1.13), p. 125], there is in fact equality in ([T7|. To go on, we need 
a truncation using the sequence a m , since the quadratic variation (M, M) t is not 
a bounded random variable in general. Then (|17p (with equality as explained 
above) and (jlip imply 

sup \M (t A T am ) — B m ((M, M) t A a m )\ 

0<t<K 

sup \M(T S A T a J - B m ((M, M) Ts A a m )\ 

0<s<(M,M)k 

sup \W(s)-B m (s)\ 

0<s<a, n A(ALM) K 

< sup \W(s)-B m (s)\. 

0<s<a, n 

Hence by Theorem [TJ with m > m^(C), 

Pi sup \M(tAT a J~B m ((M,M) t Aa m )\>ai(\og,a m )im2- 2 ?\ 

{0<t<K J 

< p( sup \W(s)-B m (s)\>ai(log„a m )im2-^\ 

lo<s<a m J 

< 10(a m 2 2m ) 1 - c . 

This proves (a). 

To show (b) it is enough to consider the inequality 

sup \M(t)-B m ((M,M} t )\ 

0<t<K 

< sup \M(t)-M(tAT a J\+ sup \M(tAT am )-B m ((M,M) t Aa m )\ 

Q<t<K 0<t<K 

+ sup \B m ((M,M) t Aa m )-B m ((M,M)t)\. 

0<t<K 

From this point the proof is similar to the proof of Lemma [U (b) . □ 

Kicfer [6] proved in the Brownian case M — W that using Skorohod em- 
bedding one cannot embed a standardized RW into W with convergence rate 
better than 0(l)n~~ (log n) 2 (loglogn)*, where n is the number of points used 
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in the approximation. Since the next theorem gives a rate of convergence 
0(l)n~3logn (the number of points used is n — K2 2m ), this rate is close 
to the best we can have with a Skorohod-typc embedding. The same remark is 
valid for Theorem H] below. 

Theorem 3. Applying the same notations as in Lemma\3\and taking a sequence 
(c m ) increasing to oo arbitrary slowly, we have 

sup \M{t) — B m ({M,M) t )\ < c m m2"T a .s. (m oo). 

0<t<K 

Under the condition of Remark^ we also have 

sup \M(t)-B m ((M,M) t )\ <g(K)? (log g(K))i a.s. (K ->■ cx^) 

0<t<AT 

/or anj/ m Zarje enough, m > m3(3). 

Proof. The statements follow from Lemma |3] in a similar way as Theorem [5] 
followed from Lemma [2] □ 

We mention that when M is a continuous local martingale vanishing at 
and there is a deterministic function / on M + such that (M,AI) t = f(t) a.s., 
then it follows that M is Gaussian and has independent increments, see |T0j 
V (1.14), p. 186]. Conversely, if M is a continuous Gaussian martingale, then 
(M, M) t = f(t) a.s., see [TD1 IV (1.35), p. 133]. In this case the "twist and shrink" 
construction of Brownian motion described in Section 2 can be extended to a 
construction of M(t) (or a simulation algorithm in practice). Namely, we have 

M{t)-B m (f(t)) <0(l)m2~T a. s . ( m ->• oo). 

Here B m (t) = 2- m S m (t2 2m ) (m > 0) denotes the nested sequence of the RW 
construction described in Section 2. 

Combining the previous results one can replace {M,M) t by the discrete 
quadratic variation N m (t) when approximating M(t) by time-changed shrunken 
RWs. 

Lemma 4. Let M be a continuous local martingale vanishing at 0, let (M, M) 
be the quadratic variation, T be its quasi-inverse (OJ) and N m be the discrete 
quadratic variation defined by U3\) . Denote by B rn the sequence of shrunken 
RWs embedded into M by Lemma [H Fix K > and take a sequence a m — 
0(m~ 7 ~ e 2 2m )K with some e > 0, where a m > K V 1 for any m > 1. 
(a) Then for any C > 3/2 and m > 7715(C) we have 

P j sup \M{t A T a J - B m (N m (t A T a J| > 2aJ i (log» a m )im2~^ 
Lo<t<i<" 

< 14(a m 2 2m ) 1 - c . 

Under the tail- condition for any C > 3/2 and m > 7715(C) it follows 

that 

P j sup |M(i) - B m {N m {t))\ > 2t4(bg„ o m )im2-? 

< 14(a m 2 2m ) 1 - c + Di^m- 1 - 
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Proof. For proving (a) we use the triangle inequality 
sup | M (t A T am )-B m (N m (tAT a J\ 

0<t<K 

< sup \M(t A T am ) — B rn ((M, M) t A a m )\ 
a<t<i< 

+ sup \B m ((M,M) t Aa m )-B m (N m (tAT a J\. (18) 

0<t<K 

Since the first term on the right hand side can be estimated by Theorem [31 we 
have to consider the second term. For m > 1 introduce the abbreviation 

L a , m = 13(Ca m log, a m )^mhl- m > 12(Ca m log, a m )imh- m + 2~ 2m . 

By Theorem [5] (a), with C > 3/2 and m > 7714 (C) it follows that 

sup \B m ({M, M) t A a m ) — B m (N m (t A T J| 

0<*<if 

< sup |S m (L((M, M) t A a m )2 2m j2- 2m ) - B ro (iV m (t A T«J| + 2" ro 

0<t<K 

< sup sup |B m (fc2" 2m ) -B m (r2- 2m )\ + 2~ m 

0<k2- 2m <(AI,M) K Aa m \r~k\2- 2m <L a , m 

< sup sup \B^}(r2- 2m )\+2- m , 

0<fc2- 2m <a m 0<r2- 2m <L Q , m 

except for an event of probability < 3(a m 2 2m ) 1 ~ c , since the difference of a 

shrunken RW at two dyadic points equals the value of some shrunken RW Bm 
at a dyadic point. 

Then we can apply estimate (@| with some C > 1 for the last expression: 
pJ sup sup \B^(r2- 2m )\ 

[o<fc2- 2m <a,„ 0<r2- 2m <L a , m 

> ^2C*'log7V supVar(BW(r2- 2m ))^ \ < 2N 1 - ' , 

where N = \a m 2 2m \ \L a , m 2 2m \ and sup fcr Var (b$ (r2- 2m )j < L a<m . Choose 
here C so that 1 — C = |(1 — C). Then a simple computation shows that 
2N 1 - ' < (a„ l 2 2m ) 1 - c , also log N < 8m log, C log, a m , and 



2C'log7Vsup Var(BW(r2- 2m )) + 2"™ < a, 4 „(log» a m )i m 2-^ 

\ k,r J 

if to > to,5(C) > TO4(C). This argument and Theorem [3Ja) applied to (fT5)l give 
(a). 

Statement (b) again follows from (a) in a similar way as in Lemma [2l □ 

Theorem 4. Using the same notations as in Lemma^ and taking a sequence 
(c m ) increasing to 00 arbitrary slowly, we have 

sup \M{t) - B m {N m {t))\ < c m m2~^ a.s. (m -> 00). 

0<t<K 
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Under the condition of Remark^ we also have 

sup \M(t) - B m {N m {t))\ < g(K)i(logg(K))i a.s. (K -> oo) 

0<t<K 

for any m large enough, m > 7715(3). 

Proof. The statements follow from Lemma [4] in a similar way again as Theorem 
[5] followed from Lemma [5J □ 

4 Independence of the DDS BM and (M, M) 

It is important both from theoretical and practical (e.g. simulation) points 
of view that the shrunken RW B m and the corresponding discrete quadratic 
variation process N m be independent when approximating M as in Theorem [4] 
This leads to the question of independence of the DDS Brownian motion W and 
quadratic variation (M, M) in the case of a continuous local martingale M. For, 
by LemmaQ] B m depends only on W and, by ([TBI . N m is determined by (M, M) 
alone. Conversely, if the processes B m and N m are independent for any m large 
enough, then so are W and (M, M) too by Theorems [TJ and [2] It will turn out 
from the next theorem that the basic notion in this respect is the symmetry of 
the increments of M given the past. Thus we will say that a stochastic process 
M(t) it > 0) is symmetrically evolving (or has symmetric increments given the 
past) if for any positive integer n, reals < s < t± < ■ ■ ■ < t n and Borel sets of 
the line U±, . . . U n we have 

p{r| ?°} = p{r-\j*}, (19) 

where T = {M{t x ) - M(s) G U x , . . .,M(t n ) - M(s) 6 U n }, T~ is the same, but 
each Uj replaced by —Uj, and J 7 " = a(M(u),0 < u < s) is the filtration 
generated by the past of M. If M(t) has finite expectation for any t > 0, then 
this condition expresses a very strong martingale property. 

Condition (|T9"|) is clearly equivalent to the following one: for arbitrary pos- 
itive integers n, j, reals < sj < ■ ■ ■ < Si < s < ti < . . . t n and Borel-sets 
Vi , . . . , Vj , XJ\ , . . . , U n one has 

P{rn A} = p{f- n A} , (20) 

where T and T~ are defined above and A = {M (si) E V\, . . . , M(sj) € Vj}. 

Our Theorem [5] below is basically a reformulation of Ocone's Theorem A 
of [S]. There it is shown that a continuous local martingale M is conditionally 
(w.r.t. to the sigma algebra generated by (M, M)) Gaussian martingale if and 
only if it is J-invariant. Here J-invariance means that M and J^adM have 
the same law for any predictable process a with range in { — 1,1}. In fact, it is 
proved there too that J-invariance is equivalent to ii-invariance which means 
that it is enough to consider deterministic integrands of the form (t) — 
I[o,r](t) ~ I(r,oo)(*)- Moreover, Theorem B there extends the above result to 
cadlag local martingales with symmetric jumps. 

Dubins, Emery and Yor in [2] and Vostrikova and Yor in |13) gave shorter 
proofs with additional equivalent conditions in the case when M is a continuous 
martingale. In these references the equivalent condition of the independence 
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of the DDS BM and (M,M) explicitly appears. Besides, in [2J, the conjecture 
that a continuous martingale M has the same law as its Levy transform M = 
J sgn(Af) dM if and only if its DDS BM and (M, M) are independent is proved 
to be equivalent to the conjecture that the Levy transform is ergodic. Below we 
give a new, long, but elementary proof for any continuous local martingale M 
that the DDS BM and (M, M) are independent if and only if M is symmetrically 
evolving. 

Theorem 5. (a) If the Wiener process W(t) (t > 0) and the non- decreasing, 
vanishing at 0, continuous stochastic process C{t) (t > 0) are independent, 
then M(t) = W(C(t)) is a symmetrically evolving continuous local martingale 
vanishing at 0, with quadratic variation C . 

(b ) Conversely, if M is a symmetrically evolving continuous local martin- 
gale, then its DDS Brownian motion W and its quadratic variation (M, M) are 
independent processes. 

Proof. To prove (a) suppose that W and C are independent. By [TUl V (1.5), p. 
181], M(t) = W(C(t)) is a continuous local martingale. For simplicity, we will 
use only three sets in showing that M is symmetrically evolving, i.e. equation 
(|2T))) holds, the generalization being straightforward: 

P {M(si) 6 Vi,M(h) - M(s) e U u M(t 2 ) - M(s) e U 2 } 

= fp{W( Xl ) e V{\ 

[ P {W(y 2 ) - W( Vl ) £U 2 -u}P {W(y} - W(y ) £ du} 
x P{C(si) G dx u C(s) G dy ,C(h) G d yil C{t 2 ) G dy 2 } 
= fp{W( Xl ) e V} 

f P {W(y} - W(y 2 ) eU 2 ~u}P {W(y ) - W(y} G du} 

JUi 

xP{C( Sl ) G dxx,C(s) G dy ,C(h) G dyx,C(t 2 ) G dy 2 } 

= P {M( Sl ) G V u M{s) - M(t} G Ui,M{s) - M(t 2 ) G U 2 } , 

using the independence of B and C on one hand and the symmetry and inde- 
pendence of the increments of Brownian motion on the other hand. 

For proving (b) we want to show that the sequences T m (k) (k = 1,2,...) 
and M(r m (j)) — M(r m (j — 1)) (j = 1,2,...) are independent. Since N m (t) 
depends only on the number of stopping times r m (fc) < t, cf. (fT5|) . while the 
shrunken random walk B m is determined by the steps 2~' m X m (j) = M(r m (j)) — 
M(T m (j — 1)), cf. Lemma [TJ this would imply their independence and so the 
independence of W and (M, M) too by Theorems [T] and [5J For this it is enough 
to show that with arbitrary integers m > 0, n > 1, < k < n, reals t%, . . . , t n 
and = ±2~ m , . . . , S n = ±2~ m (we fix these parameters for the remaining 
part of the proof) one has 

P{AnB< k nB >k } = P{Af)B< k f)B- k } , (21) 

where A< k = C\r=i{ T m( r ) — tr}, A >k is similar, but with r = k+ 1, ...n, 
A = A< k f)A >k , B< k = f)jL x {M(r m (r)) - M(r m (r - 1)) = S r }, B >k is similar, 



16 



B. Szekely Sz T. Szabados 



but with r = k + l,...n, B = B<^ fl B > i c , and finally B~ k is the same 
as Byk, but each (5j is replaced by —Sj. For, if one can reflect all 6jS for 
k < j < n without changing the probability, then one has the same probability 
with arbitrary changed signs of SjS too, since any such change can be reduced 
to a finite sequence of reflections of the above type. Let B* be similar to B, 
but with arbitrarily changed signs of SjS. Then, as we said, (|2"TT) implies that 
P{AnB} = P{4nB'}. Since P {B} = P {B*} by Lemma [TJ the desired 
independence follows. 

We will prove (|2"Tj) in several steps. 

Step 1. In condition (|19[) one can replace s by an arbitrary stopping time a 
adapted to the filtration (-7-"°): for any Uj > (1 < j < N), 

P{F\ T°}=P{F-\ J*}, (22) 

where 

N 

F = p| {M(uj M(a) eUj}, (23) 

3=1 

and F~ is the same, but each Uj replaced by —Uj. This is somewhat similar 
to the optional stopping theorem, see [TU1 II (3.2), p. 69]. Indeed, for discrete 
valued stopping times a the statement is obvious, since then 

P{F| T«}=Y j I{o = s r }V{F\ Tl), 

where {s r } denotes the range of a, including possibly oo, and 1{S} denotes 
the indicator of the set S. For every stopping time a there exists a decreasing 
sequence of discrete valued stopping times <7i almost surely converging to a. Let 
us denote the events defined according to (|2"3"|) for Oi by F% and F~ , respectively. 
Further, denote the operators projecting L 2 (fl) onto its subspace of random 
variables measurable w.r.t. J-^.. and by Pi and P, respectively. Then 

||P {Fi | ^ } - P {F | Jj} || 2 = WPJ {Fi} PI {F} \\ 2 

< ||P, (I {Fi} — I {F}) \\ 2 + ||P t I {F} PI {F} \\ 2 

< ||I {Fi} — I {F} \\ 2 + l|E(I {F} | J* f ) E(I {F} | J-°)|| 2 , 

which goes to as i — > oo. Here we used that E(I{F} | J 7 ®.) is a bounded, 
reversed-time martingale converging to E(I{F} | J-%). Hence for any e > 0, 

||P{F| ^}-P{F-\ J*}\\ 2 < ||P {F | J^} - P {Fi | Jjj || 2 

+ \\P{F-\F° a }-P{Fr\^}\\ 2 <e, 

if i is large enough. This shows that the left hand side of the inequality is zero, 
so d22]) holds. 

Step 2. Then for arbitrary reals < Uj < Vj and Borel sets Uj (1 < j < N) 
we have 

P{G| J-°} = P{G- | F° a }, 

where 

iV 

G = p) {M( Vj + a) - M(uj +a)eUj}, (24) 

3=1 
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and G is the same, but each Uj is replaced by —Uj. For simplicity we prove 
this only for two factors, the general case being similar: 

P {M(vi + <t) - M (ui + &) G U u M(v 2 + <t) - M(u 2 + a) G U 2 | 
= J I{xi ~x 2 e Ui} l{x 3 -x 4 e U 2 } 

x P{M(vi +a) - M(a) G dxt,M{ui + a) - M(a) G dx 2: 
M(v 2 + a) - M(a) G dx 3 , M(u 2 + a) - M(tr) G dx 4 | J 7 "} 

= J I{xi - x 2 G Ui} I{x 3 ~x A e U 2 } 

x P {M(a) - M(vi + a) G dxi, M(<r) - M(u x + <r) € d£ 2 , 
M(er) - M(u 2 + <r) G dx 3 ,M(a) - M(u 2 + a) G dx 4 | 
= P {M( Ul + a) - M(« x + a) G U u M(u 2 +a)- M(v 2 + a) e U 2 \ J""} . 

Step 3. Let Ar m (i) = r m (i) — r m (i — 1) and a G [0, oo). Consider the event 

S{a) = {Ar m (fc + 1) > a} 

= {inf {u > : \M {u + T m {k)) - M (r m (k))\ > 2~ m } > a} 

= ( sup {\M(u + T m (k)) - M(r m (fc))|} < 2- m l 

Lo<«<a J 

= f| {|A/(u + r m (fc))-M(T m (fe))| <2-™}. (25) 

0<w<a 

Introduce the set of dyadic numbers D\ = {r2~ l : r G Z} (Z > 0), and the events 
5j,i(o)= f| {|M( g + r m (fc))-Af(r m (fc))| < 2' m - 2^} . (26) 

gGDi ,0<g<a 

when j > m. For j fixed, {Sj l i) l _Q * s a decreasing sequence of events. Take 

oo 
1=0 

Since Sj^(a) is increasing with growing j, so is Sj(a). Put 

oo oo oo 

S*(a)= (J Si(a)= (J fl 

We want to show that S*(a) = S(a), where the latter is defined by (|25l) . 

First fix an w G 5(a). (We suppress cj in the notations below.) Then with 
this uj, 

sup {\M(u + T rn (k)) - M(r m (k))\ < 2- m } =: a < 2~ m . 

0<u<a 

If j > m is such that 2~ 3 < 2~ m — s, then us G Sj t i(a) for any Z > 0. So w G 5j(a) 
for each j large enough, consequently, uj G 5* (a) . 
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Second, fix an ui £ S(a). Then there exists a real uq < a (depending on ui) so 
that \M(uq + T m (k)) — M(r m (fe))| = 2~ m . Since the path of M is a continuous 
function, for any j > m there exists an I > and q £ D; , < q < a such that 
\M(q + T m (k)) - M(r m (fc))| > 2- m - 2~ j . That is, u <£ Sj(a) if j > m, thus 
u<£S*(a). 

In other words, we have proved above that 

{Ar m (k + 1) > a} = S(a) = lim lim Sj t i(a) 

= lim lim f| {|Af(g + r m (fc))-M(r m (fc))| <2- m -2^'}. 

gG-D i: 0<9<a 

Consequently, any event {Ar m (/c + 1) > a} = S(a) can be written in terms of 
monotonic sequences of intersections of finitely many events of the form 

{\M(q + r m (k)) - M(r m (fe))| < c} (c > 0). 

Moreover, such an approximation can be applied to {Ar m (k + 1) G (a, b]} = 
S(a) \ S(b) as well, with any < a < b. 

Step 4- First, Steps 2 and 3 imply that for any a > 0, 

P {G n S jtl (a) | J? m(fc) } = P {G- n %/(a) I T° Tm{k) } , 

because of the absolute value in definition (121)1) of Sj t i(a). Throughout Step 4 
G and G~ are defined according to (f2"4")) with <t = r m (fc), but otherwise with 
arbitrary parameters, possibly different from case to case. Then taking limit as 
j ^ oo and I — >• oo it follows from Steps 2 and 3 that 

P {G n {Ar m (fc + 1) > a k+1 } | J^ m(fe) } 

= P {G- n {Ar m (fc + 1) > a fe+1 } | J? m (fe)} • 

We want to extend this symmetry property by induction over i = fc + 
1, . . . , n + 1. Taking arbitrary reals > and integers Z > 0, r,; > 0, 
(fc + l<i<n+l) define the events 

i 

Bt= P| {M(r ro (p))-M(r m (p-l)) = «5 p }> 
P =fe+i 

i?i = P| {AT m (p)>a p }, 

p=fc+i 

i 

Xu(r) = p| {Ar m (p) S ((r p - 1)2"*, r^] } 

L. M (r) = {Si (M( n 2- 1 + ■■■ + r k+1 2- 1 + r m (k)) 

-M{n-i2- 1 + ■■■ + r k+l 2- 1 + r m {k))) > 0} , 

and B^~ , L~ t (r) similarly, but multiplying each S p by (—1). Suppose that we 
have already proved that 

P {g n Bi-. x n Hi | T° Tm(k) } = P [g~ n bt_ y n ^ | j* m(k) } , 
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where B k — fl. Define the following event as a generalization of ([25)1 : 
SijA^r) = f] {\M{q + ri 2- 1 + ■■■ +r k+1 2~ l +r m {k)) 

qeDi,0<q<a 

-Min-tf- 1 + r k+1 2~ l + T m (k))\ < 2- m - 2-i} , 
where j > m. Then by the induction hypothesis we get that 



P \ G n Bi_! n Hi n \J K^r) n L M (r) n S itjil (a i+1 ,r) \ T° m{k) 

r k+1 ,...,r i = l 



2^+1 

(fc) 



p <Gn n if,- n (J ^./(r) n L7,(r) n S idil (a i+1 ,r) I j£ 



where we agree that when r p = 2 21 + 1, the interval ((r p — 1)2~ ( , r p 2~ l ~\ in the 
definition of K^i{r) is replaced by ((r p — 1)2 ' , oo] = (2 l , oo]. Notice here that 
the events in Ki^i can be written in terms of a difference of events appearing 
in Hi, while the events in L^/ and Sijj are both of the type appearing in G, 
though Si.jj is not affected by reflections because of the absolute values in its 
definition. Then taking limit as j — > oo and I — > oo it follows that 

p {g n b % n H i+X | j* m{k) )=p{c-nB i -n H i+l \ J* m{k) } . 

This completes the induction. 

Comparing the notations introduced in this step with the ones introduced 
above, observe that B >k = B n and H >k = H n = H n+ i, if a n +i = 0. Thus one 
obtains that 

P [B >k n H >k I T° Tm(k) } = P {B- k n H >k I F° Tm{k) ) . 

Step 5. The result of Step 4 implies that 

p{A >k nB >k \T° Tm(k) } 

= J l{r m (k) + x k+1 < t k+1 ,. . . ,T m (k) + x k+ i H h x n < t n } 

xP |S >& n {Ar m (fc + 1) G dx fe+ i, . . . , Ar m (n) G dx n } | J\ m{k) ] 

= J I{r m (k) + x k+ i < t k+ i, . . .,T m (k) + x k+ i H \- x n < t n } 

xP n {Ar m (k + 1) G dx fe+ i, . . ., Ar m (n) G da;„} | ^ m(fe) } 

= p{A >4 nB;JJ? sW }. 
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Step 6. Finally, it follows from Step 5 that 

P{Ac\B} = P{A< k nA >k nB< k nB >k } 

= E(E(i{A< fc n%}i{i >t nB >fc }|^ (t) )) 

= E(i{A< k r\B< k }p{A >k nB >k I J^ m(fc) )} 

= E(i{A< k nB< k }p[A >k nB- k \F° Tm{k) )} 
= P {A< k n A >k n B< k n B- k } . 

This proves (f2"Tj) , and so completes the proof of the theorem. □ 
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